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Recent terahertz conductivity measurements observed low-power-law frequency dependence of 
optical conduction within the Mott gap of the Kagome lattice spin-liquid candidate Herbertsmithite. 
We investigate mechanisms for this observed sub-gap conductivity for two possible scenarios in which 
the ground-state is described by: 1) a U(l) Dirac spin-liquid with emergent fermionic spinons or 2) 
a nearly critical Z2 spin-liquid in the vicinity of a continuous quantum phase transition to magnetic 
order. We identify new mechanisms for optical-absorption via magneto-elastic effects and spin- 
orbit coupling. In addition, for the Dirac spin-liquid scenario, we establish an explicit microscopic 
origin for previously proposed absorption mechanisms based on slave-particle effective field theory 
descriptions. 



I. INTRODUCTION 

The spin- 1/2 antiferromagnetic Heisenberg model on 
the Kagome lattice is thought to have a spin-liquid 
ground-state due to the interplay of geometric frustration 
and large quantum fluctuations of spin-1/2 moments 1-6 . 
This model is (approximately) realized in the material 
Herbertsmithite (ZnCu3(OH) 6 Cl 2 ), which consists of lay- 
ers of Cu 2+ ions with antiferromagnetically interacting 
spin-1/2 magnetic moments arranged in a structurally 
perfect Kagome lattice 7 . There is compelling evidence 
that Herbertsmithite has a spin-liquid ground-state. It 
shows no signs of magnetic order down to the lowest ac- 
cessible temperatures ~ 30mK, well below the magnetic 
exchange scale J « 190K. Instead, at low-temperatures, 
it exhibits power law temperature dependence of spin- 
susceptibility Xs and heat capacity C v 8 , as well as a 
Knight shift that tends to a constant at low T. 9,10 These 
features are suggestive of gaplcss (or at least nearly gap- 
less) spinon excitations. Moreover, recent neutron mea- 
surements show a broad spectrum of spin-excitations con- 
sistent with a low-energy spinon continuum 11 . 

While Herbertsmithite is a promising spin-liquid can- 
didate, the detailed nature of its putative spin-liquid 
ground-state remains mysterious. Evidence for gapless 
spinon excitations led to early suggestions that the ma- 
terial might realize a U(l) Dirac spin-liquid (DSL) with 
fermionic spinons whose low-energy dispersion consists 
of two Dirac cones 12 ' 14 . Projected wave-function stud- 
ies showed that spin-wave-functions obtained from the 
DSL mean-field ansatz have very good energy with no 
variational parameters 12 (even lower energies can be 
achieved by applying a few Lanczos steps to the DSL 
wave-function 13 ). However, recent advances in density- 
matrix renormalization group (DMRG) techniques have 
enabled simulation of wide ID strips, which suggest that 
the ground state of the Kagome Heisenberg model is in- 
stead a Z 2 spin- liquid with a sizeable spin-gap 6 . The 
issue is complicated by the observation that the Z 2 spin- 
liquid found in DMRG for is very sensitive to small per- 
turbations away from the pure Heisenberg model 15 . In 
light of this observation, it is quite plausible that addi- 
tional ingredients such as disorder or spin-orbit coupling 



play an important role in determining the ground-state 
properties of Herbertsmithite. 

These theoretical difficulties highlight the need for new 
experimental probes that can more directly measure the 
properties of emergent spinon excitations. For example, 
it was recently proposed that spin-orbit coupling could 
enable neutron scattering measurements to detect the 
gapless emergent gauge fluctuations present in the U(l) 
Dirac spin- liquid theory 17 . 

AC electrical conductivity measurements offer another 
route for measuring the structure of low-energy spin ex- 
citations. Despite the host material being a strong Mott 
insulator, it was previously pointed out for a U(l) spin- 
liquid that virtual charge-fluctuations can enable power- 
law absorption deep within the Mott gap 16,19 . Indeed, 
in-gap optical absorption has been observed in organic 
spin-liquid candidates 18 , but the data does not go to 
very low frequencies. These suggestions have led to re- 
cent measurements of the electromagnetic properties of 
single-crystal samples of Herbertsmithite, at terahertz 
frequencies 20 . This work finds a low-power-law frequency 
dependence to AC conductivity inside the Mott gap, 
whose magnitude increases as temperature is lowered. 

Inspired by this work, we revisit the issue of sub-gap 
AC electrical conductivity of candidate spin-liquid states 
for Herbertsmithite. We consider two scenarios in which 
Herbertsmithtc forms cither 1) a U(l) Dirac spin- liquid 
with fermionic spinons and an emergent U(l) gauge field 
or 2) a nearly gapless Z 2 spin-liquid in the vicinity of 
a quantum phase transition to non-collinear antiferro- 
magnetic (AFM) order. Within these scenarios, we find 
three distinct mechanisms for power-law optical absorp- 
tion. All of these mechanisms produce conductivity with 
frequency dependence: a <~ lj 2 , but are expected to have 
different magnitudes. 

The first two mechanisms are specific to a U(l) spin- 
liquid and rely on linear coupling between the applied 
AC electrical field and the emergent gauge electric field. 
First, in a U(l) spin-liquid, an external electric field cre- 
ates virtual charge fluctuations that produce emergent 
electric fields. This mechanism is a strong-coupling gen- 
eralization of the Ioffe-Larkin mechanism for optical ab- 
sorption previously discussed in Refs. 16 and 19. Here, we 
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adapt strong-coupling expansion results by Bulacvskii et 
al. 21 . This approach allows us to establish a clear micro- 
scopic origin of this effect, and enables semi-quantitative 
estimates of its magnitude. 

Second, an external electric field can create a lattice 
distortion that is non-uniform within the unit cell, but 
uniform across different unit-cells. Such lattice distor- 
tions perturb the spin-system with the same symmetry 
as the virtual polarization mechanism 21 . The magnitude 
of conductivity from the first mechanism is expected to 
be parametrically smaller than that of the second, since 
virtual charge excitations are suppressed by a factor of 
(^) 3 in a strong Mott insulator (where t is the hop- 
ping and U is the on-site repulsion). However, it turns 
out that the charge-fluctuation mechanism is enhanced 
by a large numerical prefactor, and may actually be of 
the same order as the second magneto-elastic mechanism. 
Both of these effects give roughly the same order of mag- 
nitude as observed in Terahertz measurements by Pilon 
et al. 20 . 

In addition, we have identified a third mechanism for 
absorption based on magneto-elastic and spin-orbit cou- 
plings. Here, an applied electric field induces lattice dis- 
tortions, perturbing the magnetic system. Due to the 
special pattern of Dzyaloshinskii-Moriya interactions for 
the Kagome lattice, these spin-perturbations induce uni- 
form spin-currents. Absorption from this mechanism fol- 
lows a <~ w 2 ct s , where <r s is the spin- conductivity of 
the magnetic system, thereby allowing an all electrical 
probe of spin-transport. Moreover, this mechanism is 
quite generic; it is not particular to U(l) spin-liquids, 



but will also produce a 



absorption in a nearly 



gapless Z2 spin-liquid, or even in a thermal paramag- 
netic with diffusive spin-transport. Absorption through 
this spin-orbit based mechanism is suppressed by a fac- 
tor of (y) <~ 1CT 2 compared to the second mecha- 
nism described above, where D is the magnitude of the 
Dzyaloshinskii-Moriya interactions, and J is the mag- 
netic exchange coupling. 



II. OVERVIEW OF 17(1) DIRAC AND Z 2 
SPIN-LIQUID SCENARIOS 



Before describing possible sub-gap absorption mecha- 
nisms, we briefly review the slave-particle effective the- 
ory approaches for describing U(l) and Z 2 spin-liquids. 
The ground-state properties of strong (U ^> t) Mott in- 
sulators at half-filling are well described by eliminating 
doubly occupied sites in a t/U expansion to obtain an 
effective spin model i?Hciscnbcr g = JJ2(ij) s i ' Sj, where i 
and j label sites of the Kagome lattice. To describe quan- 
tum disordered spin-liquid states of this strongly coupled 
spin-theory, it is convenient to re-write the spin-degrees 
of freedom in terms of auxiliary Schwinger fermions or 
bosons: S, = \f} ta (Tabfi,b or Si = §&J i0 <r 6&»,& respec- 



Such descriptions have a U(l) gauge redundancy as- 
sociated with the local U(l) transformation — > 
(fi,bi), and naturally lead to low-energy effective the- 
ories with emergent U(l) gauge fields. 



A. U(l) Dirac Spin Liquid 

Fermionic spin-liquid states can be obtained from such 
spin-models by re-writing spins as Schwinger fermions: 

= \n a^abfi.b, subject to the constraint that the num- 
ber of fermions, n/ = 1. A spin-wavefunction can be ob- 
tained by first decomposing i?Heisenborg into a mean-field 
Hamiltonian: 



(1) 



which respects the constraint rif = 1 only on average. 
A physical spin-wave function can be obtained from the 
mean-field Hamiltonian by projecting out the unphysical 
doubly occupied sites: |* spin ) = V g \^h mf ), where V G 
denotes a projection that removes double occupancies. 
This spin- wave function can then be used as a variational 
ansatz. 

At the mean-field level, X ^ F = (fbfj.*), and A$ F = 
(fj.lfi.t)- For the Kagome lattice Heisenberg model, the 
procedure outlined above favors a state with A,j = 
and with a sign structure of \ij such that each hexagon 
has an odd number of negative bonds, and each triangle 
has an even number 12 . One can go beyond a mean-field 
treatment by incorporating fluctuations in the phase of 
Xij : Xij ~^ Xij F e 4aiJ , and in the magnitude of fi: fj, — > \i— 
a , where a and are the space- and time- components 
of the vector potential for the emergent U(l) gauge field. 

Linearizing the low-energy spinon dispersion near the 
Fermi-energy, the resulting low-energy effective theory 
consists of four-flavors of two-component Dirac fermions 
(one for each combination of valley and spin) , coupled to 
an emergent U(l) gauge field 12,14 . In imaginary time the 
effective Lagrangian density for the DSL is: 

£ = h, u [-«*> + v d t ■ k] / k)W + (jq.n)^ a£ n (2) 

Where we have introduced the 8-component field, f a , 
where a is a super-index labeling spin, valley, and Dirac- 
components, Here, r are Pauli-matrices in the Dirac 
component basis 14 , vd ~ ^|Xij F | is the Dirac velocity, 

j = ^ is the spinon-number-current with n = 

Ek.Jk-v-"/M and j qj o = "n Ek.w/M,w-!i T /k ; u, 



and 



where a and a arc continuum versions 



tively, subject to the constraint that n/ = 1 (nj = 1). of the long-wavelength components of the lattice- vector 
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potentials a and ay. 

B. Z2 Spin-Liquids 

While the Dime Spin-Liquid wave-function gives a 
good energy for the spin- 1/2 Kagome Heisenberg model 
without any variational parameters, DMRG studies in- 
dicate that a Z 2 spin-liquid with a sizeable (« 0.15 J) 
gap to spin excitations has even lower variational energy. 
For Herbertsmithite, a large-gap Z 2 spin-liquid is incon- 
sistent with observations of gapless behavior 8-10 down 
to energy scales well below 0.15 J. However, additional 
ingredients such as disorder and spin-orbit coupling, not 
present in the pure Heisenberg model, could suppress the 
Z 2 gap, leading to a nearly gapless state consistent with 
experiments. For example, slave particle mean-field stud- 
ies of large N analogs of the spin- 1/2 Kagome Heisenberg 
model suggest that the spin-1/2 model is naturally close 
to such an 0(4) quantum critical point between a Z 2 
spin-liquid and \/3 x \/3 magnetic order 22 . 

Therefore, while the main focus of this paper will be 
on the U(l) Dirac spin-liquid, we will also explore optical 
conductivity mechanisms for a nearly gapless Z 2 spin- 
liquid. In this section, we give a brief review of the slave- 
particle formulation of Z 2 spin-liquids, and their possible 
phase transitions to magnetically ordered states. 

Z 2 spin-liquids may be described in the Schwingcr 
fcrmion effective field theory approach by introducing 
non-vanishing pairing terms Ay into the spinon-mean- 
field ansatz. An alternative description can be ob- 
tained by using bosonic rather than fermionic spinons: 
= b\ a ^rbi^ subject to the constraint that ri\, = 1. 
This description also has a local U(l) redundancy lead- 
ing to emergent, compact U(l) gauge degrees of freedom. 
In the Schwinger boson formulation, the Z 2 spin-liquid 
state is a superfluid of pairs of bosons with (bi a ) = 
but (b ia bj a ) 7^ 0. The pair condensate breaks the U(l) 
gauge degree of freedom down to a residual Z 2 degree 
of freedom associated with bi — > —bi. When treated ex- 
actly, the fermionic and bosonic formulations of the Z 2 
spin-liquids are of course equivalent; and in certain cases 
one can explicitly construct a duality between the two 
descriptions, where the boson b is just a vortex in the 
fermionic paired superfluid 24 . 

Unlike the fermionic decomposition, the bosonic for- 
mulation allows a simple description of magnetically or- 
dered states as superfluid states with (&i, a ) 7^ 0. Since 
the frustrated Kagome geometry precludes collinear an- 
tiferromagnetic (AFM) order, it is natural to consider 
non-collinear (but co-planar) AFM order with (S(r)) = 
ni cos Q • r + n 2 sin Q • r, where ni • n 2 = 0, and Q is 
a lattice-scale wave-vector whose precise form depends 
on the particular type of non-collinear AFM order 23 . 
This formulation also naturally describes a continuous 
phase transition between Z 2 spin-liquid and non-collinear 
AFM, in which the spinon pairing order parameter- 
amplitude vanishes continuously without proliferating 



topological defects of the superfluid order. The critical 
theory of such a transition contains gapless, deconfined 
spinons, and the critical properties are those of a rela- 
tivists 0(4)* quantum critical point 23 . Here * denotes 
the fractionalized critical point in which the fractional- 
ized objects obey the ordinary 0(4) critical scaling, but 
in which the physical degrees of freedom are composite 
operators of these fractionalized fields. Interestingly, the 
quantum critical properties largely independent of the 
particular details of the type of Z 2 spin-liquid or non- 
collinear AFM order. 



C. Optical Absorption - Some General Symmetry 
Considerations 

On symmetry grounds, a uniform external electric 
field, E, will couple to any time-reversal even and spa- 
tially uniform operator that transforms like a vector rep- 
resentation of the Kagome symmetry group. In the ab- 
sence of spin-orbit coupling, the electric field couples 
only to spin-singlet operators. In the U(l) Dirac the- 
ory, the operators with lowest scaling dimension are 
the 15 fermion-bilinear mass terms 14 : = /r 3 /i*cr/, 
Nb = /t 3 <t/, and Nq = fr 3 fi l f, where [i % are Pauli 
matrices in the valley basis. Since Na and Nc break 
translational invariance they cannot couple to a uniform 
E, and since N B is a scalar, only its spatial derivative can 
couple to E (which would not contribute in the optical 
conductivity regime: q w 0, ui 7^ 0). 

The next lowest dimensional operators are: the emer- 
gent gauge fields e and b, and conserved spinon-numbcr 
current j f and spinon spin-current jg k (here i refers to 
the direction of the current of the k component of spin) . 
In the absence of spin-orbit coupling, only e has the right 
symmetries to couple to E. In the presence of spin-orbit 
coupling, j s may also couple to E. Therefore, we expect 
the electric field along the z-direction to couple to the 
operator Cu(i) = ^ l€% + ^ l 2 l Js> + • • • where Ai and A 2 
are, for the moment, unknown coefficients and . . . repre- 
sents operators with higher numbers of time-derivatives 
that make negligible contributions to AC conductivity at 
low frequency. As explained in more detail below, the 
scaling dimensions of e and jg are fixed by conservation 
laws, indicating that the conductivity from coupling to 
these operators scales like a <~ lo 2 . 

Similar symmetry arguments can be applied to identify 
the low-energy operators of the critical Z 2 spin- liquid the- 
ory that may couple to an external electrical field. While 
there is no gapless emergent electric field in this scenario, 
in the presence of spin-orbit coupling, the conserved spin- 
current may couple to E again giving a <~ us 2 . The emer- 
gent 0(4) symmetry of the Z 2 -AFM quantum critical 
point gives rise to 9 low-dimensional operators with the 
small scaling dimensions 23 . If coupling to these operators 
were allowed, it would contribute very low power-law fre- 
quency dependence a ~ w ' 5 based on their known scal- 
ing dimensions. However, only higher time-derivatives of 
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these operators can couple to E, giving higher power con- 
ductivity a ~ oj 2 - 5 . Therefore, the operator that couples 
to E in the Z 2 -AFM critical point is: O z . 2 = A 3 j s + . . . , 
where A3 is an unknown coefficient and . . . represent less 
relevant operators. 

These general considerations greatly constrain the op- 
erator form of the physical current operator, and predict 



conductivity scaling like a 



In the following sec- 



tions, we explore microscopic origins of these symmetry 
allowed couplings, in order to estimate the size of the 
coefficients Ai.2,3. 



III. PURELY ELECTRONIC MECHANISM FOR 
SUB-GAP ABSORPTION IN U(l) SPIN-LIQUIDS 

In states with emergent U(l) gauge fields, linear cou- 
pling between external electromagnetic fields and the 
emergent electromagnetic fields are generally allowed by 
symmetry. Dissipation by coupling to gapless emergent 
degrees of freedom was originally discussed by Ioffe and 
Larkin 16 within the slave-boson approach, who argued 
that the physical conductivity a was equal to the inverse 
sum of the boson and spinon conductivities, and a f. 



<?b + CFf 



(3) 



This result was later adapted specifically to the spinon- 
Fermi surface state and the Dirac spin-liquid by Ng and 
Lee 19 . 

Here we will revisit this absorption mechanism for the 
case of a Dirac spin-liquid in a strong Mott insulator. We 
adapt the strong-coupling expansion results of Ref. 21 to 
reveal the microscopic origin to the Ioffe-Larkin absorp- 
tion mechanism, enabling us to make semi-quantitative 
estimates of the magnitude of this effect for Herbert- 
smithite. While our results in this section are qualita- 
tively similar to the previous study of Ref. 19, we find 
some differences. In particular, we find a different t/U 
dependence of the optical absorption than obtained from 
a slave-rotor mean-field description. Moreover Ref. 19 
found that the Ioffe-Larkin mechanism produced conduc- 
tivity scaling like oil ~ u) 2 ~^ where j3 was an unknown 
anomalous dimension. Here, we argue that the emergent 
gauge invarianes exactly fixes (3 = 0. 

While real charge transitions are forbidden in the Mott 
insulating phase, virtual charge transitions can mediate 
a coupling between an externally applied physical elec- 
tromagnetic field, and the emergent gapless gauge 
and spinon degrees of freedom of the spin-liquid, thereby 
generating dissipation. In p. 3 of Ref. 21, Bulaevskii et 
al. derived an expression for the electron polarization 

and current operators for a 3-site Hubbard model trian- 

3 

gle: 1A2, in the large U limit. The resulting operators 
are expressed purely in terms of spin- variables, using a 
strong coupling expansion in t/U <C 1, to eliminate con- 
figurations with double occupancies. The leading order 
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FIG. 1. Panels a) and b) show the spin- and magneto-elastic 
mechanisms that couple the emergent gauge electric field e 
to an applied external field E shown here along the vertical 
(y) direction. Panel a) depicts a spin configuration that has a 
net polarization along the electric field due to virtual charge 
fluctuations. Dark shaded bonds labeled + and light shaded 
bonds labeled — respectively indicate strong and weak spin- 
singlet correlations. Panel b) shows a distortion of the Cu 2+ 
ions induced by the electric field. This distortion alters the 
bond distance, which increases the exchange coupling Jij for 
dark bonds labeled + and decreases it for light bonds labeled 
— . Both mechanisms have the same symmetry. 

contribution to the polarization for a triangle of spins 
was found to be: 



S 3 -(S 2 -Si) 

r«„ . 
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where e is the charge of an electron (which should not 
be confused with e, the emergent electric field), a is the 
lattice spacing (which should not be confused with the 
emergent vector potential a). Here, V is the discrete 
lattice derivative defined on the lattice of bond-centers 
with bonds labeled by their endpoints i and j. The x 
and y components of this discrete gradient are written 
out explicitly in the first line. 

This polarization will couple the spins of each triangle 
linearly to an external electric field E, giving a perturba- 
tion to the spin-Hamiltonian: 



SH E (t) 



J d 2 rP(r,t) -E( 



r,t) 



(5) 



A small magnitude applied AC electric field E = 
Eoe _iLJt produces an induced polarization: (P)e — 
(P u P- u )E(ui), delivering energy through SHE at a rate 
of 3?e [(P U P_ W )] Eq, or, equivalently, giving conductivity 
a — iw(P w P_ w )riA where ha = ^ s the density of 

triangles in the kagome lattice with lattice spacing a. 

Expressing the polarization operator P as a sum of 
the polarizations on each triangle, re-written in terms 
of spin-operators through Eq. 4, we find that the spin- 
electric field coupling produces the conductivity: 



o-ilM 



7f w(12e)2 



u 



<|v(s, • s,)H 



(6) 
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In the large U limit, the bond energy gradient is di- 
rectly related to the emergent electric field (see Ap- 
pendix A for derivation): 



J 



V(Si • s,) 



(7) 



where the second line is Faraday's law of induction for 
the emergent gauge field, which becomes asymptotically 
exact at low energies due to the emergent gauge symme- 
try. In a deconfined U(l) spin- liquid state with emergent 
scale invariance, such as the Dirac spin-liquid, the electric 
field cannot acquire any anomalous dimension due to this 
emergent conservation law. These considerations exactly 



where J w is the magnetic exchange coupling. From nx the sca li ng f 0y l ~ ^ 2 for the Dirac spin-liquid. 



this expression, we see that there is a linear coupling be- 
tween the physical electrical field and the emergent elec- 
tromagnetic field. Such a linear-coupling was previously 
identified by Ioffe and Larkin based on a slave-boson ef- 
fective field theory description of the Hubbard model. 

The result presented here is simply a strong-coupling 
generalization of the Ioffe-Larkin result valid for t <C U. 
Though the functional form of the E • e coupling is identi- 
cal to the slave-rotor treatment, the strong-coupling ex- 
pansion predicts a different t/U dependence compared 
to a mean-field slave-rotor treatment of the Hubbard 
model. For example, comparing to Eq. 6, and using 
a = ^j(jphiph), we identify the physical electromagnetic 
current (in the optical conductivity limit, q = 0, w ^ 0) 



as: 



J P h 



(8) 



This result disagrees with the slave-rotor mean-field the- 
ory in the large U limit, which would have predicted 
j P h siavc-rotor ~ {p>d t e (see Appendix B for more details). 
We believe that this discrepancy is due to the mean-field 
treatment of the slave-rotor theory, which does not prop- 
erly account for the gauge constraint in determining the 
properties of virtual charge excitations. A simple illus- 
tration of the shortcoming of the slave-rotor mean-field, 
is that it predicts that gapped charge excitations propa- 
gate with velocity ta ~ vp rather than J a as would be 
appropriate for the strong-Mott limit (see Appendix B 
for more details). In contrast, the strong-coupling t/U 
expansion maintains the physical electron Hilbert space, 
and in our view, provides a more reliable determination 
of the coupling between spins and external electromag- 
netic fields. The results of the t/U expansion can then 
be taken as a starting point for the slave-particle effec- 
tive field theory, which we can better expect to correctly 
capture the low-energy physics of the spin-system. 

These considerations show that the contribution to AC 
conductivity from coupling to the emergent gauge field 
is: 



72tt e 2 t 2 , 



(9) 



In the absence of magnetic instantons the emergent gauge 
field components form a conserved current: 



7 M = 

J cm 



b 

ex z 
d t b + z ■ V x e = 



(10) 



To estimate the magnitude of the prefactor, we need 
to rely on an approximate treatment of the low-energy 
effective field theory. In a random-phase approximation 
(RPA) treatment of the low-energy effective action, the 
gauge field propagator is set by the spinon electromag- 
netic response kernel. Starting with the low-energy ef- 
fective Dirac theory, Eq. 2, choosing a gauge such that 
a = 0, and adding a source term, J e , for e ku) = — u;ak, w : 



£ =/k,« + v d t ■ k] / k)U , + j q ,n • a q ,n 



■ wa ka 



(11) 



Integrating out the spinons within the RPA yields the 
following effective action for a: 

£ a = \ai^K)\^ to)a\_ u - ■ .a ku (12) 

Where Kj — (j'p) is the spinon electromagnetic re- 
sponse kernel. Integrating out the emergent gauge field 
gives the effective action for the source J e : 



2cr D ' 



I J, 



k,wl 



(13) 



where an is the conductivity for massless Dirac fermions 
(with the effective charge set to 1), which is a universal 
constant at low-frequency. Analytically continuing back 
to real time, these manipulations show that, in the RPA, 
the real-time (ee) correlator takes the form: 



i i i \ t0J c 

(e„eiJ R p A = -^T^j 



(14) 



The precise value of this universal constant can de- 
pend on whether one has free, interaction dominated, 
or disorder-dominated fermions. Neglecting the effects 
of gauge- fluctuations, one has: a D = ^ = |. 25 
Combining this result with Eq. 6 yields: 



ctil(w) w 48^^ 



IF J h 



(15) 



Here the large numerical prefactor is primarily due to 
the factor of 12 in the relation between spin and electric 
polarization, Eq. 4. 

The experiments by Pilon et al. 20 , examine optical 
conductivity at terahertz frequencies corresponding to 
huj w 4meV« 50K. At to = ITHz, they measure 



6 



three-dimensional conductivity of order 1730 (ITHz) ss 
0.1f2 _1 cm _1 . For interlayer spacing ~ IOA, this cor- 
responds to two-dimensional conductivity er(lTHz) » 

10 -4 ^. Comparatively, we have U ~ leV and t ~ U/10. 
These numbers give an estimated magnitude for the con- 
ductivity produced by coupling to the emergent electric 
field: oil (ITHz) ~ 10~ 5 ^-, which is slightly smaller 
than, but roughly the same order as that observed. We 
will see in the next section, that a coupling with the 
same symmetry also arises from magneto-elastic cou- 
pling, which is expected to contribute additively to this 
purely electronic contribution. 



IV. MAGNETO-ELASTIC MECHANISMS FOR 
ABSORPTION 

In the previous section, we discussed sub-gap optical 
absorption in a U(l) spin-liquid through the Ioffe-Larkin 
coupling of the applied and emergent electric fields. Be- 
yond the pure spin-model, there will be additional path- 
ways for optical conduction from magneto-elastic cou- 
pling between spin and lattice distortions and from the 
interplay of magneto-elastic and spin-orbit couplings. In 
practice, we will argue that these mechanisms may be 
more important than the Ioffe-Larkin type contributions. 
Furthermore, in the presence of spin-orbit coupling, these 
additional absorption mechanisms will also be relevant 
for nearly gapless Z2 spin-liquids without an emergent 
U(l) gauge field. 




FIG. 2. The Kagome point group is generated by mirror 
reflections, M x , about the x-axis passing through the center 
of the vertical hour-glasses, and n/3 rotations, R n /3 about 
the center of the hexagons. 

Here, we represent the resulting pattern of JJ^ 's picto- 
rially for a single unit cell consisting of one upward and 
one downward facing triangles. Black (white) circled- 
bond representing SJij > (6Jij < 0) respectively, and 
it is implicit that the pattern shown should be repeated 
uniformly throughout the Kagome lattice. 

As illustrated in Fig. 2, the Kagome point-group is 
generated by 1) clock- wise | rotations about the cen- 
ter of a Hexagon, denoted R„/3, and 2) mirror reflec- 
tions through the horizontal axis (x-axis) passing through 
the points of triangles, denoted M x . The transformation 
properties of the basis elements in Eq. 17 are: 



M,. 



Magnetoelastic Coupling to Emergent Electric 
Field 



Since the Cu ions have a net charge, an external electric 
field can give rise to a relative displacement of Cu ions 
within the unit cell. As shown in Fig. lb, this distortion 
will increase the hopping strength along some bonds and 
decrease it along others, yielding a perturbation: 



SH } 



ME 



(16) 



While deriving the precise distortion resulting from an 
electric field would require knowledge of complicated de- 
tails of various lattice-elasticity parameters, it is sufficient 
for our purposes to simply determine what spin-operators 
can couple to this distortion on symmetry grounds. The 
pattern of modified exchange couplings, SJij, induced by 
the lattice-distortion shown in Fig. lb can be expressed 
as linear combinations of the basis elements: 



= S4 • S5 — Si • S2 
= S4 • S3 — S3 • S2 
= S5 • S3 — S3 • Si 



(17) 



This three dimensional representation of the Kagome 
point-group reduces into a one-dimensional pseudoscalar, 
-B2, irreducible representation (irrep) spanned by: 



Or 



1 

71 



(19) 



and a two-dimensional vector, E\, irrep spanned by: 



o\ 



o% 



1 

71 
1 

71 



(20) 



The external electric field E is a vector, transforming 
like the E\ irrep of the Kagome point-group, and hence 
only couples to the E\ components: SHme ~ E • Oe x - 
Comparing the symmetry of Oe x , we see that this opera- 
tor is just the discrete bond-energy gradient correspond- 
ing to the polarization operator generated described in 
the Eq. 4 of the previous section. Therefore, as previously 
pointed out in Ref. 21, the magneto-elastic coupling due 
to the electric field induced displacement of Cu ions will 
generate a coupling of the same form as in Eqs. 4,6, but 
with a coupling constant f^-^gf ~ K ^ \i in place of 

jj) 3 appearing in Eq. 4, where Uij is the ionic dis- 



the 
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placement induced by the electric field, and Kq u is the 
effective spring constant for the Cu displacement shown 
in Fig. lb. 

For the Dirac spin-liquid, the perturbation Eq. 16 pro- 
duces a linear coupling between the external and emer- 
gent electric fields, just as with the Ioffe-Larkin mech- 
anism discussed in Sec. III. In this case, following the 
derivation of Eq. 15 in Sec. Ill, we see that this magneto- 
elastic coupling leads to optical conductivity: 

In contrast, since this magneto-elastic mechanism relies 
on the presence of an emergent U(l) gauge field, it is 
not expected to produce low-power-law frequency depen- 
dence of optical absorption for the nearly critical Z 2 spin- 
liquid scenario. The energy scale Kc u a 2 is expected to 
be of order ~ leV, which is comparable to U. Conse- 
quently, for the Dirac spin-liquid, this absorption mecha- 
nism may dominate over the Ioffe-Larkin contribution of 
Eq. 15, which is smaller by a factor of (^) . 

B. Spin-Orbit Coupling 

The mechanisms described in the previous section re- 
lied on the existence of an emergent electric field, a low- 
dimensional operator with the right symmetries to couple 
to the external electric field (i.e. a time-reversal invari- 
ant, spatially uniform, spin-singlet operator transforming 
under the Ex representation of the Kagome point group) . 

These mechanisms are special to the U(l) spin-liquid, 
since, in the alternate scenario where the ground state 
is a nearly critical Z 2 spin-liquid, there are no spin- 
singlet operators with the right symmetry properties. 
However, as we will show in this section, the presence 
of spin-orbit coupling enables the external electric field 
to couple not only to spin-singlet operators, but also 
time- reversal invariant spin-pseudo- vectors. In particu- 
lar, the combination of spin-orbit and magneto-elastic 
couplings will generically enable the electric field to cou- 
ple to spin-current, enabling purely electrical probes of 
spin-transport. This coupling to spin-current provides al- 
ternative mechanisms for optical absorption which lead 
to AC electrical conductivity scaling like a <~ uj 2 in any 
gapless relativistic state. In particular this mechanism 
is also relevant in a nearly gapless state, at finite fre- 
quency above the spin-gap, for example in the vicinity of 
an 0(4) quantum critical point between a Z 2 spin- liquid 
and non-collinear antiferromagnet. 

1. Dzyaloshinskii-Moriya (DM) Terms 

We begin with a brief overview of spin-orbit coupling 
in Herbertsmithite and related Cu spin- 1/2 Kagome an- 
tiferromagnets. Due to the atomic spin orbit coupling in 



the Cu d-orbitals, A so , and because spin-exchange occurs 
along a bent Cu-O-Cu super-exchange pathways, Her- 
bertsmithite has substantial Dzyaloshinskii-Moriya (DM) 
interactions: $^<»j) ^ij ' (^» x Sj), where the DM vector 
D w sin a where a is the angle of the Cu-O-Cu oxy- 
gen bond (defined such that a = for a linear super- 
exchange pathway). For Herbertsmithite, the DM inter- 
actions are expected to be roughly 10% of the exchange 
energy: D rts 0.1J. 

The dominant effects of the DM interactions come from 
the z-component of the DM vectors, z ■ Dy, where the z 
denotes the c-axis perpendicular to the Kagome planes. 
As we will see below, the DM terms generate coupling 
between external electric fields and spin-current. Since 
the in-plane components sum to zero within the unit cell, 
these components can only create a response to spatially 
inhomogeneous electric fields that vary in sign over short 
distances of order a single lattice spacing. Such sharp 
spatial variations arc not relevant for AC conductivity 
measurements that probe the response to long wave- 
length fields, and therefore the in-plane DM components 
may be neglected for what follows. 

In a Schwinger-fermion mean-field description of the 
Heisenberg model, the DM terms simply become spin- 
orbit hopping terms for the spinons: 

/ 

D.; ' (Si x S,) -+ i «Dy • (4<W jl6 ) (22) 

where — (flafj^)- For the Dirac spin-liquid state, 
this pattern of DM interactions translates into a quan- 
tum spin-Hall mass for the Dirac spinons. Such a term is 
expected to gap out the Dirac spinons, enabling gauge- 
instantons to proliferate. Due to the quantum spin-Hall 
response, such instanton flux-insertion events are bound 
to spin-flips, and their proliferation leads to magnetic 
order 14 . Similarly, Schwingcr-boson mean-field treat- 
ments of the Kagome-Heisenberg model indicate that 
the Z 2 spin-liquid state for the bare Heisenberg model 
is destabilized by finite DM interactions, which are pre- 
dicted to drive the model to non-collinear antifcrromag- 
netic order 22 . Experimentally, Herbertsmithite appears 
to avoid magnetic order to the lowest achievable temper- 
atures, and is either gapless or has a very small gap. 



2. DM Modulation 

The lattice distortion discussed above in Sec. IV A al- 
ters the Cu to Cu bond lengths following the pattern 
shown in Eq. 20, thereby changing the hopping strengths 
along these bonds in the same manner. Since the DM vec- 
tor on bond ij is related to the Cu-Cu hopping strength 
by: Dij ~ , this distortion also changes the magni- 
tude of the DM vectors in the same manner, introducing 
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a) 



b) 




FIG. 3. Blue and green dots are Cu and O ions respectively, 
a) DM arises from the bent Cu-O-Cu bonds. The pattern of 
D vectors is shown by ® and indicating D into and out 
of the plane respectively for the bond-orientation shown with 
arrows, b) In an electric field, the oppositely charged Cu 2+ 
and O 2- ions undergo relative displacements. The O ion dis- 
placement changes the bond-angle, and the Cu displacement 
changes the bond exchange strength as described in Fig. 1. 
As explained in the text, this alters the DM pattern as indi- 
cated by the + and — labels, which couples the electric field 
to the spin-current. 



the perturbation: 



SH 



Cu Distortion 
DM 



5^ ay* • (S< x (23) 

(ij) 



For definiteness, we consider E along the vertical (y-) 
direction as shown in Fig. 3b. In this case, the magni- 
tudes, |ay|, are twice as large on the horizontal bonds 
(12 and 45) as on the angled bonds (23, 31, 34 and 53). 
The sign structure of ay can be represented pictorially 
by: 



DM 



2X - X. - 



(24) 



where we have expanded the coupling as a linear combi- 
nation of the basis operators: 

+ 

%=z- (Si x S 2 ) + z- (S 4 x S B ) 

+ 

+ ^ + =z- (S 3 x S 2 ) + z- (S 4 x S 3 ) 

+ ^ + =z- (S 5 x S 3 ) + z- (S 3 x Sx) (25) 

which form a closed set under the action of the Kagome 
point group. 

Since the oxygen and Cu ions are oppositely charged, 
they will generally they will undergo a relative displace- 
ment in a uniform electric field, E. In particular, the 
oxygen bond angles will be changed thereby changing the 
DM terms as indicated in Fig. 3b. To linear order the 
DM change is affected only by the bond-angle change, 
which depends only to the component oxygen displace- 
ment perpendicular to the local bond direction within the 
kagome plane. For definiteness, we consider E along the 



vertical direction in Fig. 3 resolving the applied field in 
to directions along and perpendicular to the local bonds: 
-E/y,|i and Sy,x, as shown in Fig. 3b. It can be readily 
seen that this contributes a term of the same form as 
Eq. 24. 

This perturbation, 5H® M , can couple to any spin- 
operator with the same symmetries. Therefore, it is in- 
structive to determine its transformation properties of 
under the Kagome lattice point-group by considering ir- 
reducible representations constructed from the basis of 
Eq. 25. In the presence of spin-orbit coupling terms, the 
Kagome point-group remains unchanged but the spatial 
point-group transformations must also be accompanied 
by spin-rotations 14 . Since, as explained above, we may 
neglect the in-plane components of the DM vectors the 
only modification is that mirror reflections about the x- 
axis must be accompanied by 7r spin-rotations around the 
x-axis. 

Under the action of and M x , the vector chirality 
combinations of Eq. 25 transform like: 



+ 

{ X ' X+ ' +X 1 



M,. 



+ 

-{Y, + YvY + } (26) 



where the minus sign in the first line arises because 
rotations invert our chosen bond-orientation, swapping 
Si x Sj — > Sj x Si — —Si x Sj. Similarly the minus sign 
in the second comes from the extra spin-rotation that ac- 
companies the mirror reflection M x . Furthermore, since 
all of these combinations arise from linear combinations 
of vector spin-chirality, each of the basis vectors are in- 
variant under time-reversal symmetry. 

The three dimensional representation of the Kagome 
point-group spanned by the basis shown in Eq. 25 is re- 
ducible, and decomposes into one-dimensional pseudo- 
scalar, A2, irreducible representation (irrep) spanned by 
(see Appendix C of Ref. 14 for more details about the 
representation theory for the Kagome point-group): 



O 



1 



DM.Ao 



12 

and a two-dimensional vector, E\, irrep spanned by: 



(27) 



O v — — - 
~ v/6 



n x 



V2 



(28) 



The Cu bond-, and oxygen bond angle- distortions cou- 
ple only to the E x vector piece: 5H DM ~ E • O v . To 
summarize, Yltij) a ij^ ' (Si x Sj) is time-reversal invari- 
ant, transforms as the z-component of a vector under 
spin-rotations, and transforms like a E± vector represen- 
tation of the Kagome point-group. These transformation 
properties are exactly the same as the z-component of the 
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spin-current, jgz . Therefore, there one should generically 
expect a linear coupling between the electric field and 
the spin-current, mediated by the oxygen bond distor- 
tion shown in Fig. 3b. Interestingly, this coupling implies 
that one can directly probe spin-conductivity of Kagomc 
materials with the DM pattern shown in Fig. 3a, purely 
through electrical measurements. 

This expectation can be explicitly verified. Noting that 
the spin-current operator obeys the continuity equation 
Q§r- = —i[H,S a ] = — V • J 5a, and using the Heisenberg 
Hamiltonian to compute the time-evolution, one finds 
that the a-component of spin-current along the ij bond 
is just (Sj x Si) a . Alternatively, for the spin-liquid sce- 
narios, one can perform a mean-field decomposition of 
the perturbation term within the Schwinger fermion (or 
boson) framework using Eq. 22. In either case, one finds 
that the perturbation, SH VM ~ E O^, induces a uniform 
spin current, jgz, perpendicular to the applied E field in 
the z x E direction. Hence, the optical conductivity re- 
sulting from this mechanism will be proportional to uj 2 
times the spin-conductivity. 

Modeling the Cu and Oxygen bond-angle distortions 
as springs with effective spring-constant Kq u and Kq re- 
spectively, one can estimate the magnitude of the induced 
coupling: 

W D „«-^E.(ixj Sl ) (29) 

JiVeffd 

where K cS = ( Kq \ 0& - — j^r) is tnc effective spring 

constant accounting for the fact that the Cu 2+ and 2 ~ 
ions have opposite charge, and a is the Cu-0 bond-angle. 

For the Dirac spin-liquid, Eq. 29 indicates that an ex- 
ternal electric field induces a spin-current js* ■ 

]s* ~ ttt (iuxtd) (z x E) (30) 
J K e sa 

where <jd is the spinon-conductivity. Similar results will 
hold for the Z2 case. 

Since jg* is a conserved quantity, its scaling dimen- 
sion is fixed in any scale invariant theory, like the Dirac 
spin-liquid or in the vicinity of a quantum-critical point, 
resulting in <j{uj) ~ oj 2 . 

aDM ra i (7) (£^) as (31) 

where o s = j^{js'js') is the spin-conductivity. 

For the Dirac spin-liquid: a s w or) is constant. 
Comparing to the contribution from the spin-symmetric 
magneto-elastic coupling of Eq. 21, we see that the mag- 
nitude of the Ioffc-Larkin mechanism may dominate due 

to the factor of (y) w 10~ 2 . However, for the scenario 
of a nearly gapless Z 2 spin-liquid, the spin-conductivity 
coupling may be the dominant absorption mechanism. 



3. Spin-Orbit Mechanism in a Thermal Paramagnet 

While we have illustrated this result explicitly for the 
U(l) Dirac spin- liquid, the DM mediated coupling be- 
tween electric field and spin-current is generically al- 
lowed by symmetry, and will be present in any theory 
of Kagome spin systems with DM interactions. Further- 
more, for any scale invariant theory, the resulting con- 
ductivity from coupling to this phonon mode will scale 
as ~ w 2 , including a Dirac spin- liquid or any nearly gap- 
less Z-i spin-liquid in the quantum critical regime near a 
continuous phase transition, for example to non-collinear 
antiferromagnetic order. 

Furthermore, the coupling of an electric field to spin- 
current may be the most important mechanism for ther- 
mal paramagnetic states at intermediate to high temper- 
ature ranges, well above the low-energy scale Ag. For 
example, the low energy scale Ag would be the spin-gap 
of a Zi spin-liquid, the magnetic ordering scale of the 
frustrated antiferromagnet, or the spinon bandwidth in 
a 17(1) spin-liquid (note that this scale is well-separated 
from J in a frustrated lattice like the Kagome). Regard- 
less of the nature of the spin ground-state at T <C Ag, 
for T > Ag one can expect diffusive spin-transport with 
non- vanishing spin-diffusion coefficient D s , which tends 
to D s w J a 2 at high temperature 26 . D s is related to 
spin-conductivity, a 8 , by the Einstein relation a s « D s x s 
where Xs is the spin-susceptibility. At intermediate 
to high temperature one expects Curie- Weiss-like spin- 
susceptibility: x s ((x),T) rts t+Ocvj w ith #cw ~ J > 0. In 
Herbertsmithite, for example, Xs is a slowly increasing 
function of T down to T w 50K 8 , giving rise to a constant 
spin- conductivity, and power-law electrical conductivity 
a <~ lo 2 . We note that this effect may be difficult to ob- 
serve in practice, even at higher temperatures, due to the 
suppression of this conductivity-pathway by a factor of 

V. DISCUSSION 

To summarize, we have found two different mecha- 
nisms for sub-gap absorption in a spin-liquid Mott in- 
sulator on the Kagome lattice. For a U(l) spin-liquid, 
absorption can occur through coupling the external elec- 
tric field to the emergent electric field mediated either 
by virtual charge fluctuations or magneto-elastic effects. 
This mechanism produces a <~ uj 2 conductivity for the 
Dirac spin-liquid. 

Second, the interplay of spin-orbit and magneto-elastic 
couplings creates a direct coupling of an electric field to 
the spin-current of the Kagome system. Therefore, op- 
tical conductivity measurements are also probing spin- 
conductivity. This mechanism is quite generic, and will 
produce a ~ to 2 electrical conductivity in any gapless 
relativistic phase or critical point, including the Dirac 
spin-liquid and the vicinity of the 0(4) deconfined crit- 
ical point between a Z2 spin-liquid and a non-collinear 
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magnetically ordered state. 

In contrast, this mechanism will not lead to substan- 
tial power-law conductivity deep inside a magnetically 
ordered state. A magnetically ordered state can be re- 
garded as a spin-superfluid, with spin-waves correspond- 
ing to the superfluid density mode. In this state, the 
spin-current is given by the spatial gradient of the super- 
fluid phase, and due to the huge mismatch in the velocity 
of light, c, and the spin-wave velocity, Ja, an electromag- 
netic field can only couple to multi-spin-wave emission 
processes giving a negligible contribution. 

Therefore, at asymptotically low frequency and tem- 
perature, power-law optical conductivity can likely be 
taken as a sign of gapless (or nearly gaplcss) spinon ex- 
citations. At intermediate to high temperatures where 
the spin-system is a thermal paramagnet, the coupling to 
spin-current will produce a <~ uj 2 due to spin-diffusion, 
regardless of the nature of the low-temperature spin- 
phase (its magnitude may be small due to the factor of 
(y) in Eq. 31). However, at temperatures below ~ 50K 
the spin-susceptibility \ s dips sharply with decreasing 
T, 10 indicating the onset of correlations and a departure 
from the trivial thermal paramagnet. 

Within the two spin-liquid scenarios, the temperature 
dependence of the coefficient of the <~ uj 2 optical ab- 
sorption is more subtle. Within the scaling regime for 
both theories, the optical conductivity will act like co 2 
times a universal function of ^. For the nearly critical 
Z2 spin-liquid, the dominant mechanism is proportional 
to the spin-conductivity. For the 0(4)* critical point, in 
the high-frequency {uj » T) regime relevant to the ex- 
periments of Pilon et al. 20 , the conductivity is expected 
to gradually increase, saturating to a constant value as 
^ — >• 00 (see e.g. Ref. 27). This temperature dependence 
is in agreement with that observed in Ref. 20. 

As described above, in the Dirac spin-liquid case, the 
Ioffc-Larkin type mechanisms dominate. The result- 
ing conductivity is proportional to the spinon-resistivity 

Pd(^) = • Here the finite temperature and 

frequency scaling is known only in the limit of large- 
number of Dirac flavors, Nf. For Nf — > 00, one recovers 
the case of free Dirac fermions for which po decreases 
with decreasing temperature in the uj 3> T limit. This 
Nf = 00 prediction shows the opposite temperature de- 
pendence of that observed in Pilon et al.'s experiment 20 . 
However we caution that it is unclear to what extent this 
prediction applies to the physical case of Nf — 2, given 
the absence of a controlled theory in this limit. For ex- 
ample, in other strongly coupled field theories based on 
holographic duality the temperature dependence at high 
frequency can be of arbitrary sign depending on details 28 . 
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Appendix A: Microscopic Identification of the 
Emergent Electromagnetic Fields 

The emergent gauge degrees of freedom arising in the 
slave-rotor effective theory are sometimes described as 
"fictitious" , as they arise out of a redundancy of the slave- 
rotor formalism. However, while the emergent vector po- 
tential, a M , has no direct physical meaning, gauge invari- 
ant quantities derived from a M have observable physical 
meaning in terms of spin. 

The functional form of the expression of the emergent 
electric and magnetic fields may be deduced on symme- 
try grounds and for a strong Mott insulator. Further- 
more, their parametric dependence on the parameters 
of the Hamiltonian can be fixed by dimensional analy- 
sis, at least for a strong Mott insulator where there is 
only one length scale, a, and one energy scale, J. It is 
nevertheless gratifying to produce an explicit derivation 
of the microscopic meaning of the gauge flux through 
an elementary plaquette and electric field along a lat- 
tice link. One should keep in mind, however, that these 
microscopic operators will be renormalized in going to a 
low-energy effective description and will receive contribu- 
tions from all other low-energy operators with the same 
symmetries. 

By considering spatial Wilson loops of a M within the 
slave-particle effective theory, Refs. 29 and 30 demon- 
strated that the emergent magnetic flux through a tri- 
angular plaquette corresponds to scalar spin-chirality of 
the triangle: 

bi jk ~ S, • (Sj x S fc ) (Al) 

where sites k form corners of the triangle (travers- 
ing cyclically in the clock- wise direction). As shown in 
Refs. 21 and 32, in a strong Mott insulator, the vec- 
tor spin-chirality is just the electrical current circulating 
around the triangle in the clock-wise sense: 

bijk ~ m « e ) S * ' (Sj x S fc ) + O (j^j 

(A2) 

In this appendix, we follow a similar route to establish 
the physical meaning of the emergent electrical field e. 
We will show that, in terms of the physical spin-degrees of 
freedom, the electric field corresponds to the gradient of 
bond-energy, which in turn corresponds to the electrical 
polarization arising from virtual charge fluctuations (as 
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described in Ref. 21). 



the identities: 



Spatial Wilson Loops and the Emergent 
Magnetic Field 



We begin by reviewing the arguments of Ref. 29 
to identify b as the scalar spin-chirality. In the 
Schwinger-fermion approach one rewrites spins as aux- 
iliary fermions: Sj = j} fi,b, subject to the con- 
straint J2a fiafi,*? — 1- Converting the resulting Heisen- 
berg Hamiltonian into a functional integral with fields 
/, / (corresponding to /t , /) and introducing Hubbard- 
Stratonovich fields Xij to decouple the resulting four- 
fcrmion hopping term yields the following term in the 
effective action: 

£j =J fi,a( T )fjA T )fj,b( T )fhb{ T ) 

^ JXij fi,afj,a + J\.ij fj.afi.a + l^XijXij (A3) 

In the resulting low-energy theory, fluctuations of the 
magnitude of Xij , around their saddle-point (mean-field) 
value are massive, and the phase fluctuations about the 
saddle point value are simply e laij . The Wilson loop 
around a spatial loop V is: 



Wr 



IF 

or 



or 



(A4) 



where a,j is the vector-potential on the bond (ij). 



This Wilson loop can be generated by introducing an 
infinitesimal source term 

SjC = fjijXij + VijXij (A5) 
By making the change of variables: 

Xij -> Xij - Vij Xij -> Xij - fjij (A6) 

one shifts the source for Xij to a source for f} a fj, a , indi- 
cating that the Wilson loop is equivalently expressed in 
terms of a string of fermion fields: 



Hi II'"' ^\[f,.„f,„ 

or or 



(A7) 



For example,_the flux through a triangular plaquette is 

^iAj fi,afj,afj,bfk,bfk,cfi,c- 

One can then restore the Grassman expression to oper- 
ator form 31 : W^, -> fi ia fj,afl b fk,bft t Ji,c which can be 
expressed purely in terms of physical variables of through 



ft c i o i < - 

Ji,aJifi = —Tlf,i + »» ' (Tab ~ — + Oj • (T a b 



2 2 

fi,afi.b = "H 2 ~ n f,i) - S * ' a o,b 



Sab 



i,b 2 v_ "J i*/ ~' - <*u ■ - 2 



Si • (Tab 

(A8) 



where is the spinon occupation number of cite i, 
which can be regarded as approximately constant: n/ w 
1, and (j are the spin- 1/2 Pauli matrices. 

The magnetic flux through a triangle ijk is then = 
i (W ijk - W^) = 3m (Wy fc ) « S, • (Sj x S fc ). The 
same procedure identifies the magnetic flux through a 
square plaquette with vertices ijkl as: 



%fe( — g (^y'fc + ^jfe/ + bkli + hij) (A9) 



2. Emergent Electric Field 



Proceeding by analogy to the previous section, the elec- 
tric field on a link (ij) is given by the space-time Wilson 
loop: 



1 



= lim — 

J st^o St 



a i(oy(t)+oJ(t)-ay(t+(5t)-aJ(t)) 



S, • ^Sj x Sj^ Sj ■ ^Si x 



+ h.c. 



(A10) 



where we have explicitly forced to be Hermitian in 
order to avoid operator ordering ambiguities in restoring 
the path-integral expression to operator form. The time- 
evolution of Si follows from evolving under the Heisen- 
berg Hamiltonian H = JJ2{ij) S» ■ Sj: 



Si = -i[H, S] = J^S J+d xS, 



(AH) 



where d are the vectors connecting site i to its nearest 
neighbors. Combining Eqs. All and A10 gives: 



Si • Sj+dj - Sj ■ s i+dt 



di 



(A12) 



where d i} dj arc the nearest neighbor vectors of sites i 
and j respectively. 

For a 3-site triangle, Eq. A12 reduces to a discrete 
gradient of bond-energy (on the lattice of bonds), which 
was shown by Bulaevskii et al. to be proportional to the 
polarization arising from virtual charge fluctuations. 
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(i,t +St) (j ,t+5t) 

5 ^ / 

t* > \ 

(i,t) a,t) 
— ^ r ij 

FIG. 4. Space- Time Wilson loop corresponding to the emer- 
gent electric field dj on lattice-bond (ij) . 

3. Physical Meaning of the Emergent Gauge Field 
in the Presence of DM Interactions, and Field 
Induced Ordered States 

In the presence of DM interactions, the above deriva- 
tion of the physical meaning of b and e is altered. As 
described previously, since the in-plane DM vectors van- 
ish within the unit cell, the uniform z— component of 
the DM interactions is expected to be the most impor- 
tant. On general symmetry grounds, the DM terms allow 
b ijk -+ S t ■ (Sj xS k )+a (Sf + SJ + Sg), where a is of 

order: a»0(f). The precise form of a may be derived 
from the Wilson loop derivation outlined above but incor- 
porating the DM interactions at the slave-fermion mean- 
field level as a spin-orbit coupling term for the spinons. 

The appearance of this new term shows directly that 
any perturbation that couples to the z-magnetization will 
tend to induce an emergent-gauge-magnetic flux. For ex- 
ample, as discussed in further detail in Ref. 29, these 
considerations demonstrate that neutron scattering ex- 
periments can in principle be used to detect gapless fluc- 
tuations in b. 

Furthermore, it has implications for the behavior of 
Herbertsmithitc in an applied magnetic field. For ex- 
ample, the Zeeman term ^f-H z from an external field 
H z along z will induce a magnetization M z = XsH z , 
which will produce a corresponding gauge magnetic flux 
bind ~ &xsH z . In perfectly clean Dirac spin- liquid, the 
density of states vanishes at the Dirac point, and the spin 
susceptibility grows linearly with H: Xs' rac ~ tt- Em- 
pirically, however, Herbertsmithite appears to have spin- 
susceptibility that tends to a constant of order a fraction 
of Xs ~ j a t low temperature 8-10 . Within the U(l) 
spin-liquid scenario, one could possibly interpret this fi- 
nite susceptibility as a finite density of states induced, 
for example, by impurity scattering. 

In this case, the induced magnetization per spin from 
an external field H z would be: M z ~ producing an 
emergent gauge flux through each triangle of order $t = 
bo 2 ~ ^ - 10~ 3 H'T' 1 . This emergent gauge flux 

has a magnetic length £ b ~ 30a^/ ~ 10 2 ^\/r^E' 
which corresponds roughly to the ordinary orbital mag- 
netic length an electron gas in the field H z . Therefore we 
expect that in a U(l) spin-liquid, the spinons will see this 
induce 6 flux as an orbital field, with an effective 6-flux 
density that is comparable the physical flux density of 



H z . This is effect is surprisingly large, since charge mo- 
tion is largely frozen deep in the Mott insulating phase, 
and purely electronic mechanisms give effective flux for 
spinons that is smaller factor of (^) <~ 10~ 3 compared 
to the physical flux 32 . 

This large orbital b field, would then induce Landau 
levels into the Dirac spectrum. Neglecting Zeeman split- 
ting, and in the absence of interactions would have a four- 
fold degenerate zeroth Landau level, with an approximate 
SU(4) symmetry as in single- layer graphene 33,34 . How- 
ever, the infinite density of states will generically cause 
interactions to remove this degeneracy and split the ze- 
roth Landau level, gapping out the spinons and allowing 
instantons to proliferate. The nature of the resulting bro- 
ken symmetry state will depend on the manner in which 
the Landau-level symmetry is broken, as different break- 
ings of the SU(4) symmetry result in different quantum 
Hall responses (e.g. quantum- spin, valley, etc... hall 
effects). Within the Dirac spin-liquid scenario, such ef- 
fects could possibly account for the field induced spin- 
freezing past a critical external field strength of a cou- 
ple Tesla, which has been observed in Herbertsmithite 
through NMR measurements 35 . Moreover, this scenario 
could also be relevant in a weakly gapped Z 2 spin-liquid 
scenario, where the Zi state arises out of the Dirac spin- 
liquid state by BCS pairing. In such a scenario, the min- 
imal field required to induce spin-ordering could possibly 
be interpreted as H c for the spinon superconductor (one 
would naively expect spinon superconductors formed by 
adding pairing to a U(l) spin- liquid to be strongly Type- 
1 due to the large effective gauge-charge for the emergent 
gauge field). 

At strong-magnetic fields where the Zccman-cncrgy de- 
termines splitting of the lowest Dirac Landau-level, there 
will be a quantum spin-Hall response 36 , tying a spin- 
flip to each instanton insertion leading to a field-induced 
magnetically ordered state. At lower fields, alternative 
scenarios might occur 34 . For example, if the zeroth 
Landau-level splits by valley polarization, there will be a 
resulting quantum valley-Hall effect. This would endow 
instantons with sub-lattice symmetry breaking quantum 
numbers; and proliferating such instantons would result 
in a valence-bond crystal (VBC) state. This raises the 
possibility of a various sequences of continuous phase 
transitions as a function of field, for instance from Dirac- 
spin-liquid to low-field VBC state, to a high-field mag- 
netically ordered state. 

The emergent gauge electric field will also be altered by 
the presence of DM interactions. On symmetry grounds, 
e can receive contributions of order a from operators like 
(Sj x Sj), which correspond to spin-current in the spinon 
theory. However, unlike the magnetic field case described 
above, it is difficult to externally couple to spin currents, 
so it is unclear whether these extra terms enable useful 
new probes. 
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Appendix B: Slave Rotor Effective Field Theory 

In this paper, we have used the strong-coupling t/U 
expansion of Bulacvskii et al. 21 to incorporate virtual 
charge fluctuations into the effective spin-model descrip- 
tion of the Hubbard model. An alternative approach that 
was previously used in Refs. 16 and 19, is to work directly 
with the Hubbard model using a slave-boson or slave- 
rotor effective theory. However, as we will show here, 
while qualitatively similar, the mean-field treatment of 
the slave-rotor theory does not capture the correct t/U 
dependence of charge-excitations in a strong Mott in- 
sulator. To show this fact, we begin by reviewing the 
slave-rotor framework and its predictions for the electro- 
magnetic response of a U(l) spin-liquid Mott insulating 
state. 

The slave-rotor approach begins by re-writing the 
physical electron operator at site i and with spin a as 
a product of a spin-less boson bi and spinful fermionic 
spinon f ia : c i-IT = bifi ta , subject to the constraint: 
rib = n/ — 1. It is technically convenient to describe 
the boson as a U(l) quantum rotor: bi — e l0i . This 
description suffers from a U(l) gauge redundancy cor- 
responding to the local transformations: Oi — > 0, + A j; 



fi, 



-*A. 



fi a . As a result the corresponding effec- 



tive theory will contain a fluctuating compact U(l) gauge 
field, of. 

By a choice of gauge, one can associate the physical 
electron charge with b, and a Mott insulating state is 
then described by a charge gap A for the bosonic rotor- 
excitations. The resulting low-energy effective theory 
reads: 

S = J dr(L b + L f ) 

\{d T + a° i +A < () fti| 2 + A 2 |6i| 2 



2U 

J^t^e^+^bibj 



(ij) 

L f = J2h( d T-V- a °) h E >i/" iJ Mi ( B1 ) 

» (ij) 

where A 1 / is the physical electromagnetic field, /i is de- 
termined such that {rib) = (n/) — 1, and where we have 
treated bi as a soft-rotor by relaxing the rigid constraint 

N 2 = i. 

The hopping amplitudes t' and t b can be estimated 

tij(fifj) and 
t independent of in- 



within a mean-field approximation: t\j 



i{0i 



Whereas t 



teraction strength U, the spinon bandwidth decreases 
stronger interactions starting from « t in the vicin- 
ity of the Mott transition (A <C t,U), and reducing to 
tJ — > J w jj for a strong Mott insulator (U 3> t). In 
the large U limit, we see that the slave- rotor effective 
action simply reproduces that of the Schwinger fcrmion 
treatment of the Heisenberg spin-model, without discard- 



ing the gapped charge fluctuations that enable non-zero 
optical conductivity. 

Here we see the first hint that the mean-field treatment 
of the slave-rotor theory does not quantitatively capture 
the dynamics of charge excitations: the mean-field treat- 
ment predicts that charge excitations propagate with ef- 
fective mass 2^2, whereas, in a strong Mott insulator one 
would expect strong coupling between spin fluctuations 
to enhance the holon mass to ~ ( j^t)- 37 



a. Speed of Emergent Light 



The Maxwell term for a generated by integrating out 
the gapped boson fields in the slave-rotor action of 
Eq. BI, takes the form ee 2 + -6 2 , where e = — \7a + id T a 
and b = z ■ V x a are the emergent electric and magnetic 
fields respectively, and e and fi are the corresponding di- 
electric constant and magnetic permeability. 

We can compute the speed of emergent 'light', c* = 
-^=, within the slave-rotor theory for a strong Mott in- 
sulator. The dielectric constant e can be determined by 
examining the action cost of a uniform electric field cor- 



responding to ao = 0, a 



E 



<9fc.£fc 



-iQr . 



<9 fe >£fc 



uj.k 



h tfe 2U k 



-#?(0,0) 



3» 2 
UP 



2U 



E 



de k 
dki 



t z a' 



n 2 



(B2) 



where we have kept only the leading order contributions 
in t/U and identified A sa U, as appropriate for a strong- 
Mott insulator. In the first line, we have used the fact 
that the diamagnetic contribution to Kb is equivalent 
to explicitly subtracting if°°(0,0); this just ensures that 
Kb(0, 0) = as required by gauge invariance. This result 

f 2 2 

identifies: e <~ -jjr- 

The magnetic permeability, yu is obtained from 
K l b J '(Q = 0, q — > 0), which by gauge invariance is propor- 
tional to q 2 — q % q* . Therefore, it is sufficient to compute 
only the K™ component: 



E 



(de k y 
V9fe' J 



k L 2U + £ k+q/2\ [ 2U + £ k-q/2\ 



Kg (0,0) 
(B3) 
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Expanding the denominator to 0(q 2 ) gives: 



(ds k \ 



1 de k 3 dsk dei 

4 dk'dk™ i 4 gfcj 3fc"- 



2\2 



,fe I — 2t7 — i U 2U ) \-^2U-) 

r>f»vj [£— 1 

k J . 



-12U 4 q l q m 



(lu 2 + u 2 y 



(B4) 



where the first term in the second line is the average of 
the band-mass over the full band, which vanishes, and 
Vi = ||t is the velocity. For (nearly) isotropic bands, we 
can approximate the velocity average by ~ v F Si m , and 
we find that 



jtf(n = 0,«-K))~<{^ 



(B5) 



This gives: e 
results gives: 



-, and - 



(taf 
U 3 



Combining these 



ta <~ Vf 



(B6) 



where vf is the bare Fermi velocity corresponding to half- 
filling and hopping t. Again, in a strong Mott insulator, 



one would expect the only relevant scale to be the spin- 
exchange J, not t, further fueling our suspicions that the 
slave-rotor mean-field is not sufficient for a quantitative 
description of charge dynamics. 



b. Physical Electric Current Operator 

The physical electrical current is defined as the vari- 
ation of the effective lagrangian density with respect to 
the electromagnetic vector potential: j p h = Fol- 
lowing the previous section, the effective action for the 
emergent gauge field in the presence of an external phys- 
ical electrical field, after integrating out the gapped bo- 
son fields takes the form: Ce = e(e + E) 2 , which con- 
tains a linear coupling of the emergent and physical elec- 
tric fields: Cec — 2ee • E, which, in the gauge where 
E = —d T e is equivalent to (after integrating by parts): 
Cec — — 2e<9 T e • A. Meaning that the physical electrical 
current is: 



t z d z 

j p h = 2eed T e ~ e-^-d T e 



(B7) 



Comparing to Eq. 8, we see that the slave-rotor theory is 
a factor of t/U different than the strong-coupling expan- 
sion, which we attribute to the aforementioned shortcom- 
ings of the mean-field slave-rotor treatment of the gapped 
charge dynamics. 
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